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ABSTRACT 

X-ray observations of galaxy groups and clusters are inconsistent with the 
predictions of the simplest hierarchical clustering models, wherein non-baryonic 



Q . and baryonic components are assembled together under the sole influence of grav- 

—I. . ity. These departures are in the sense that the intergalactic medium is hotter 

and more extended than expected, and become increasingly strong for less mas- 
sive systems. I model these effects by constructing baseline sequences of hydro- 
, ^ , static poljd;ropic models normalized to observations of high-temperature clusters 

t^ ! and numerical simulations, and then transforming them by adding proscribed 

amounts of heat per particle at the cluster center. I present sequences with a 
universal value of this heating parameter that simultaneously reproduce recently 
Q>^ '. published observed (gas and total gravitational) mass-temperature and entropy- 

r-| ! temperature relations. The required amount of energy injection is consistent 

^! with constraints on the number of supernovae needed to account for observed 

2 '. intracluster silicon abundances, provided that energy injection is centrally con- 

to '. centrated. I argue that most of the heating occurred during or after the assembly 

of the cluster, and not exclusively in pre-collapse proto-cluster fragments. 



C^ I Subject headings: galaxies: clusters: general — intergalactic medium — cosmol- 

ogy 



1. Context 

Hierarchical clustering in a universe where non-baryonic dominates over baryonic matter 
is the primary heuristic framework for studying large scale structure formation; and, one of 
the primary laboratories for disconfirming variations deriving from particular combinations 
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of fluctuation spectrum and cosmological world model is the ensemble of clusters of galaxies.0 
Such distinguishing characteristics of the cluster population include the two-point correlation 
function and mass function of collapsed non-baryonic halos (e.g., |Governato et al. 1999| ). 



Since non-baryonic matter is not amenable to observation, these relationships must 
be studied indirectly through gravitational lensing and baryonic tracers such as stars (in 
galaxies) and hot intergalactic gas. The latter, perhaps the most unambiguous tracer because 
of the density-squared dependence of the hot gas emissivity, is the focus of the present study. 
The now well-developed database of cluster X-ray observations made throughout this decade 
by the i^O^^Tand ASCA observatories have provided accurate measurements of the cluster 
temperature function (e.g., pVlarkevitch 199^ ), and log N - log S relation (e.g., |Jones et al. 



1998| ); however, mapping from these to the relevant non-baryonic halo relations requires 
knowledge of the thermal and djTiamical structure and history of the intracluster medium 
(ICM). 

In the most well-studied hierarchical clustering scenarios, structure formation is merger- 
driven, proceeding sequentially from small to large scales. A remarkable result of numerical 
simulations along these lines is that the mass density profiles of relaxed, non-baryonic halos 
can be characterized by a universal functional form ( INavarro, I'renk, fc White 1997| - here- 



after NFW; IJing 1999| ) - at least over some large radial range. These halos are not precisely 
self-similar, since less massive halos tend to collapse earlier when the universe was more 
dense and thus have - on average - more concentrated non-baryonic density distributions. 
Subject to coUisional processes such as shocks, the structure of the baryonic component will 
differ from that of the non-baryonic, but in a regular and predictable manner if the two 
components simply approach equilibrium in concert (e.g., |Eke, Navarro, fc Frenk 1998| ). 

However, there is ample observational evidence for departures from such simplicity - 
evidence that becomes more striking as one focuses on less massive systems. The most 
well-established such breakdown is the relative steepness of the correlation between X-ray 
luminosity and ICM temperature (T^,), with L^ oc Ta.~^ observed in contrast to L^ oc T^."^ 
predicted if gravitational collapse is the sole mechanism driving the evolution of both com- 
ponents (e.g., Arnaud fc Evrard 1999|) . Likewise, values of T^ very near the non-baryonic 



component virial temperature are predicted, an equality that evidently increasingly breaks 
down for less massive systems based on departures from firupa'^ = kT^ (e.g., [Bird, Mushotzky" 



k Metzier l995| ) and M oc T,^/^ ( |Horner, Mushotzky, fc Scharf l999| ; [Ettorl t ^'ablan 1999 



^For simplicity - and because there is no universally accepted demarcation - I follow Jones et al. (1998) 
in referring to all bound collections of galaxies as clusters, therein subsuming the usual definition of groups 
of galaxies. 
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- hereafter EF), where ^rrip is the mean mass per particle, a the one-dimensional cluster 
velocity dispersion, and M the total mass (evaluated at v^qq: the radius within which the 
average density contrast, relative to the critical density, is 500). There is also evidence for a 
concomitant increase in ICM extent from the steepness of the Mjcm-Tx correlation (where 
MicM is the ICM gas mass within rsoo; |Mohr, Mathiesen, fc Evrard 1999| ), and the gradual 



flattening towards low temperatures of the relation between entropy and T^ (the "entropy 
floor;" Ponman, Cannon, fc Navarro 1999| - hereafter PCN). 



All of these deviations from predictions where ICM structure is exclusively determined 
by gravitational collapse are indicative of an important role for extra-gravitational heating 
in determining ICM properties and their variation with mass. Understanding the magnitude 
and source of the heating not only is necessary for mapping from observed X-ray properties 
to those of the non-baryonic component, but also sheds light on galaxy formation since 
feedback from the earliest star formation epoch is the most likely heating mechanism. Since 
most metals in the universe are produced in objects of high overdensity (|Cen fc Ostriker 



1999| ), and are subsequently ejected into intergalactic space by supernova-driven galactic 



outflows ( Pettini 1999|) , measurements of ICM heavy metal enrichment provide independent 



constraints on the magnitude of energy injection ( [White 1991| , [Loewenstein fc Mushotzky] 
IT9961). 



I address this set of interconnected puzzles in the present work as follows. Adopting the 
most widely-used functional form for the non-baryonic density distribution that is completely 
determined by a single concentration parameter, baseline mass sequences of equilibrium, 
polytropic ICM configurations are constructed. The sequences are anchored at the high 
mass end using the results of high-temperature ICM observations and simulations (that are 
mutually consistent), and then extended to lower masses following specified mass-scalings 
of the non-baryonic halo concentration. A thermodynamic perspective proves physically 
illuminating ( [Bower 1997| ; [Balogh, Babul, fc Patton 1999| ), and so heated families of mass 



sequences are calculated that differ from the baseline sequences by fixed amounts of heating at 
the cluster center. It then remains to be seen whether the mass sequence in a particular family 
can simultaneously explain all of the observed deviations from the scaling relations predicted 
by models without heating, while also conforming with the nucleosynthetic constraints from 
ICM abundance studies. I will show that, in fact, they do so. 



2. The Baseline Configuration 

In this section I construct simple analytic models that accurately reproduce the observed 
and simulated ICM distribution of the most massive clusters. This serves as a useful point 
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of departure for creating heated distributions for the entire range of cluster masses. 



2.1. Polytropic Models 

I consider spherically symmetric, hydrostatic gas configurations in a total gravitational 
potential following the NFW prescription. Thus, 

1 dPgas _ G'Mgrav(<r) 



2 



Pgas dr r 

where pgas and Pgas are the ICM density and pressure distributions, respectively. Mgrav is 
the total gravitational mass given by 

iVfgrav(< r) = Mvirial ^ TT , (2) 

where 

X 

FNFw(a;) = ln(l + x) - -— — ; (3) 

1 + X 

Anew is the non-baryonic mass distribution scale-length, c the concentration parameter 
'^viriai/ctNFW) a^d r^iriai the virial radius. Models that are analytic and conveniently flexi- 
ble are obtained under the polytropic assumption, Pgas oc Pgas^^^''" QSuto, Sasaki, fc Makino| 
19981 ; [Cavaliere, Menci, fc Tozzi 1999| ; EF). The polytropic gas distribution is determined by 



the single function 



-1 



X 



(4) 



The gas temperature and density follow from Tgas = To^(r/aNFw) and pgas = Po^'^ ('"/ctNFw)) 
respectively, where To and po are the (assumed finite) central values of gas temperature and 
density. The constant 

i"NFw(c) {n + l)To ^ ' 

must not exceed a maximum value to assure non-negative temperature and density out to 
^virial, ^- c there is a minimum ratio of To/T^m&\ for each pair {c,n). In the isothermal limit 
{n -^ oo), the density distribution is given instead by 

Pgas = Poo6'iso('^/aNFw), (6) 

where 

^iso(x) = (l+x)^-/^ (7) 
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A. 



/C -tvirial /q\ 

Fnfw(c) Tgas 

( PVIakino, Sasaki, fc Suto 1998| ; EF). The gas density at the origin for the isothermal case is 
related to that at infinity by po = Poo^^'^° ■ 

The virial mass, radius, and temperature are connected through the relationships 

-'"virial „ '"virial '-^cPcrit \y) 

and 

"^-'virial J- ^-''^ virial ( ^ r\\ 

pTflp / '"virial 

where Ac is the density contrast (with respect to the critical density pcrit = 9.21 10^'^°/i7o^ g 
cm~^; hjQ is the Hubble constant in units of 70 km s~^ Mpc~^) within which virial equilibrium 
is established. 

In addition to the concentration parameter, c, and polytropic index, ra, I parameterize 

models by quantities more observationally relevant than To and po - the global gas fraction, 

/virial = /gas ('"virial) 5 and dcusity-squared- Weighted average temperature in virial units, r. It 

follows that 

To hnic) 



rT^'^^^ (11) 



-^virial -'2n+l (c) 

and 

^ -cVvirial/n(c)-\ (12) 



Po ^c 3. 



Pcrit 3 

where 

/^(c) = f'dx x^e(xr. (13) 

Jo 
Equation (11) must be solved iteratively. 



2.2. The Baseline Model 

In order to establish a baseline sequence of models without extra-gravitational heating, 
I first turn to the literature on numerical simulations. In particular, I consider the average 
z = model in an Qo = 0.3, Aq = 0.7, hjo = 1 cosmology from Eke et al. (1998). Although 
the demographic properties of the cluster population greatly depend on world model, the 
z = structure of the most massive clusters is relatively insensitive to this choice. The 
non-baryonic component in the simulated cluster is completely characterized by Mviriai = 
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1.51 lO^^hjo^^ Mq, a, = 100 (r^iriai = 2.99/i7o"^ Mpc, fcT^Mai = 6.9 keV), and c = 6.47; the 
gas component by /viriai = 0.087 (by assumption), T^ ^ T^inni ij = 1)? and 



^2 \-3/3flt/2 



Pgas ~ Pcrit^o 1 H o ) (14) 



2 
core 



with (5o = 1970, rcore/?"viriai ~ 0.05, and ^gt = 0.735. 

Temperature profiles measured with ASCA have been characterized by n ^ 5 polytropic 
distributions ( [Markevitch et al. 1998 ); however, both temperature ([Irwin, Bregman, fc 



Evrard 1999| ) and density profiles in hydrostatic models with this index are significantly 



steeper than those in the simulations. In fact, the density distribution in the ICM of the 
simulated cluster is accurately recovered only in models approaching the isothermal limit. 
Moreover, the best fit has r = 1.1 rather than r = 1. Figure 1 shows a comparison between 
the simulated cluster density profile expressed by equation (14) and models with /viriai = 
0.087 and (ri, r) = (5, 1.1), (100, 1), and (100,1.1) - clearly the last provides the superior 
fit. The cause of the deviation of the temperature profile in the numerical simulation, with 
its mild gradient at large radii, from its characterization in the baseline model is evidently 
due to residual unthermalized bulk kinetic energy in the simulated ICM at r > 0.5rviriai- 
This has little effect on the observed cluster-wide temperature - generally, an emission- 
weighted average over ~ 0.5rviriai ^ that is the the global temperature used in the correlations 
considered in this paper. 

In a recent paper, EF describe the results of their analysis of ROSAT PSPC data on 
the ICM in 36 high-temperature clusters. The mean parameters that emerge from surface 
brightness fitting, using functions of the forms expressed in both equations (6) and (14), 
display a remarkable agreement with the above. They find mean values of /3fit and aNFw/''^core 
of 0.72 and 3.17 compared to 0.735 and 3.08 in the simulations of Eke et al. (1998). However, 
while (on average) A^^o = 11-3 in the latter, EF find Ai^o = 10.3 from fitting the surface 
brightness profile using equation (6) - corresponding to r = 1.1 as expected (Figure 1). The 
mass-temperature relationship derived by EF is also consistent with r = 1.1 at the high-Tj. 
end. 

The appropriate value of the global gas fraction, /viriai, can be estimated from the values 
of MicM at the highest T^ (where effects of heating should be minimized) calculated within 
rsoo by Mohr et al. (1999). Q Extrapolating to rviriai using an n = 100 polytrope with r = 1.1 



•^ICM masses derived in EF (1999) generally agree with those in Mohr et al. (1999) for clusters in common 
to the two papers, although there are significant discrepancies in gas mass fraction due to differences in 
adopted temperatures and mass estimation technique. 
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in an NFW potential with rviriai/aNFW = 6.2 corresponding to kT^ = 10 keV (Mvinai = 
2.3 lO^^/iro"^ Mq) yields /vMai = 0.155/i7o"^/^ (extrapolated from /ga^ = 0.136/i7o"^/^ at 

^500 )• 

To summarize, I adopt as baseline configurations hydrostatic models with a polytropic 
index of n = 100 normalized such that the emission-measure-weighted average temperature 
is 1.1 X the virial temperature, and gas fraction within the virial radius fvinai = 0.155/i7o~^'^ 
(po = 5050/i7o~^'^Pcrit)- The gravitating mass is assumed to follow equations (2) and (3), 
with 

where Mio = 2.3 lO^^h^Q^ M0 corresponds to the virial mass at a density contrast of 100 for 
a /cTj, = 10 keV cluster, and Cio = 6.2 is adopted from Eke et al. (1998) assuming 7 = 0.1 
(appropriate for a CDM-like fluctuation spectrum; NFW). These models yield extremely 
accurate representations of the density distributions of observed average high-temperature 
clusters (as well as of numerical simulations of the same), and do so for the correct observed 
mean temperature. Any proposed heating mechanism must therefore have relatively little 
effect on the ICM of the most massive clusters. 

The adopted models are virtually isothermal. Models where the inverse of the poly- 
tropic index is substantially different than zero (or where the temperature distribution is 
otherwise strongly varying with radius) have density distributions that are too steep to be 
consistent with the observations (as well as the simulations) if the ICM is in hydrostatic 
equilibrium in an NFW potential (this remains true for other concentration parameters than 
those considered above). Thus, if the temperature profiles reported in Markevitch et al. 
(1998) (recently brought into question by Irwin et al. 1999 and [Kikuchi et al. 1999]) hold up. 



either hydrostatic equilibrium, electron-ion equipartition ([Ettori fc I'abian 1997|; [I'akizawa 



19981 ), or the assumption of a NFW- type total mass profile must break down well within the 



virial radius. 



3. Heated Families of Polytropic Mass Sequences 

I now describe a simple prescription for constructing families of mass sequences for 
polytropic intracluster media in hydrostatic equilibrium in an NFW potential that deviate 
from the baseline sequence described above by having been heated by differing amounts. 
Heated families are completely determined by specifying the final gas fraction and polytropic 
index for each cluster mass. I show, for several different families, that a mass sequence 
corresponding to the heating needed to explain the observed entropy floor also reproduces 



the other departures from gravity-only self-similarity described in Section 1. 



3.1. Formalism 

Consider a transformation from the baseline polytropic configuration to a new one: 

K,T^p^'}-^{n,T„,pJ, (16) 

where the bl superscript refers to the baseline model parameters. The changes in specific 
entropy and heat corresponding to the transformation are 



As 



3 k 

as = -- 



and 



Ag 



Tds 



3 k 



2 firUp 



2 firrip 



rp rpbi 



, To 2 po 

ln--77 In^T 

T^' 3 p^' 



(17) 



3To-T^U d\nT 



dlnp 



dT 



where ds = (3k/2fimp)d\n{T / p"^^^) and the integrals extend from the baseline to the new 
configuration. One can define an effective adiabatic index, F, such that 






po_ 



r-i 



from which it follows that 



h T 
As = Cp 1-— In^ 



prup 



and 



where 



and 



k 



piTLp ^ 



rphl 



rp rphl 



Cr 



2r-2 



Dr = 1 + Cr (r - 1)-' 1 



3r-5' 
r-1 



d\np 



For a constant heating adiabat 



dlnp 1 

dXnT ~ F-r 



To-T^U dlnT 



Dt 



dT 



(19) 

(20) 
(21) 
(22) 
(23) 



(24) 
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T.\F^ 



T, = T'J + CrT,-p, = pi'[l + Crjl^) , (25) 

where Tg is defined such that, at the center, the new configuration has an extra heat per 
particle, Ag = kTg/finip relative to the baseline configuration. Non- decreasing entropy 
requires F < 1 (temperature increasing as density decreases) or F > 5/3 (temperature 
increasing as density increases). F = 1, F = 0, and F = oo correspond to isothermal, 
isobaric, and isochoric transformations, respectively; for F = 1 equation (25) is replaced by 

T, = T'J;p, = p'^e-^. (26) 

I assume that the amount of heating per particle is the same for all clusters (and therefore 
less effective for more massive clusters), so that 

rc = 7f7^ = e{——\ , (27) 

Jvirial V ^Wio / 

where e is the (assumed universal) injected heat per particle in units of 10 keV. To make a 
clearer connection with observationally relevant quantities, in place of Tg and F, e and the 
slope z/ are specified, assuming that the global gas fraction in the heated configuration is 
given by 

/virial = /lO (^^l^j , (28) 

where /lo is the global gas fraction for a cluster with a 10 keV ICM (/viriai is generally 
assumed constant with mass along baseline mass sequences). 

The prescription for calculating heated families of polytropic mass sequences can be 
summarized as follows. (1) Baseline mass sequences are determined from specification of the 
concentration-mass relation (the parameter 7 in equation 15), /lo (0.155/i7o"^'^), n*' (100), 
and r (1.1). (2) A heated family is then computed, having specified u and n, by varying e. 
In practice, for each mass (and its appropriate /viriai from equation 28 and concentration c 
from equation 15) hydrostatic models are iterated on the parameter Ar = (To — Tj''')/T^ina.i 
until the desired value of e is achieved. The relationship between e and Ar is approximated 
as 

3^Mviriai\^^^^, 2 HpJpI^) \ ^29) 



2 V Mio / V 31n(l + Ar/r)/ ' 

which is exact if the heating is along an adiabat with a constant index F. 
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3.2. A Simple Heated Family 

For simplicity, I initially consider a heated family where the gas fraction, /viriah ^ind 
polytropic index, n, are identical to those in the baseline model {i.e. z/ = and n = 100; 
also 7 = 0.1). Perhaps the most striking evidence for significant heating of the ICM in a 
form similar to that of equation (27) is the entropy floor at low T^ recently reported by PCN, 
where the "entropy" is defined as kT^nf,'"^^^; rtf. is the electron density in cm~^ evaluated at 
Tc {tc is defined as one-tenth the radius within which the total density contrast is 200, as 
estimated from T^ using the baseline sequence of models), and kT^ is in keV. The value of 
e that produces an entropy-temperature relation with the observed entropy floor at ^ 100 
keV cm^ is 0.35.0 This (dashed) curve is shown in Figure 2 along with observed data points 
(assuming hm = 1) from analysis of ASCA imaging data ( |l:''ukazawa 1997| ; |Hwang et al. 



119991 ). Also shown is the value from numerical simulations (filled triangle; Eke et al. 1998), 



and the (solid) curve for e = 0. Note that the latter is slightly steeper than linear due to the 
more concentrated non-baryonic mass distributions for cooler clusters. 

Remarkably, this same simple heated mass sequence is in excellent agreement with all 
of the other observed trends alluded to in Section 1. Figure 3 show a plot of total mass 
evaluated at t^qq versus T^ for e = (solid curve) and e = 0.35 (dashed curve). While most 
of the heating in the e = 0.35 sequence goes into decreasing the gas concentration, rather 
than increasing the temperature, less massive clusters are predicted to become increasingly 
hotter than expected in the absence of heating. For this sequence, F ^ 0.8 and Cr ~ 0.15 
(see equations 19 and 22), nearly independent of mass, so that all clusters have an increase in 
internal energy per particle of C^kTq ^ 0.5 keV. The resulting predicted mass-temperature 
relationship is 

M(r5oo) ^ 1.32 10^^ (^^^^^^V^)' ^ ^(^)^™''Me. (30) 

Q{M) is a form factor relating M(r5oo) to Myiriai and increases from 1 at kT^ = 10 keV 
to 1.12 at kTx = 1 keV for 7 = 0.1 in equation (15). Equation (30) is in excellent accord 
with estimates from X-ray observations compiled by Horner et al. (1999; data points in 
Figure 3) and the best-fit mass-temperature trend from EF (stars in Figure 3). A similar 
accordance is found for the temperature dependence of gas mass within v^qq (Figure 4, where 
the stars represent the best-fit trend from Mohr et al. 1999). In a similar vein to Figure 3, 



''This value is sufficiently large that is has a non-negligible effect on even the hotter clusters, and f have 
reduced the dimensionless average temperature r in the baseline model from f .f to f .05 to preserve the high 
level of agreement with the observed density profiles for high-Tj; clusters. 
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Figure 5 shows the e = and 0.35 model variations of {kT^ij-ii^i/ fiiripY^'^ with T^ that can be 
characterized by 

^-'virial ~ 7777 , (31) 

1.U5 

compared to observed values of the velocity dispersion a from Girardi et al. (1998; data 
points) and the best-fit a-T^ correlation of Bird et al. (1995; stars). The agreement with 
the best-fit trend is excellent; although, a large observed scatter is evident, and I have 
not considered velocity dispersion anisotropics and gradients that are possible sources of 
discrepancy between kT^i^ai and finipa'^. 

To summarize, if one considers cluster non-baryonic matter and gas distributions as 
predicted by large scale structure numerical simulations, and then transforms the gas distri- 
bution by adding 3.5 keV per particle of heat at the center while preserving the global baryon 
fraction and polytropic index, the observed trends wherein the ICM becomes increasingly 
hotter and more extended^ are simultaneously and accurately reproduced. 



3.3. Other Heated Families 

With families of mass sequences anchored at high-T^ as described in Section 2.2, there re- 
mains the freedom to vary the following three (one for the gravitating mass, two for the ICM) 
slope parameters: the mass-scalings of total matter concentration - 7 = — dlnc/dlnMviriai - 
and gas fraction within the virial radius - u = dln/vinai/dluMvinai, and the polytropic index 
in the heated configuration - n = dlnpgas/dlnTgas. The implications of varying these are as 
follows. 

Since an increase in the concentration of the gravitating mass distribution increases the 
central ICM density in the equilibrium configuration, a steeper dependence of c on Mvinai (as 
would be expected if the fluctuation spectrum is flatter than in CDM cosmogonies; NFW) 
would imply that more heating (larger e) be required to explain the observed entropy floor. 
For example ~ 4.5 keV per particle would be required if 7 were increased from 0.1 to 0.2 
(that is the concentration increases by ~ 2.3 instead of ~ 1.5 for 1 keV, relative to 10 keV, 
clusters). 

Positive values of the parameter u correspond to cooler, less massive clusters having 



^Mohr et al. (1999) claim an absence of evidence of a more extended ICM in cooler clusters because of the 
lack of a clear trend of the average gas particle location with temperature; however, the scenario described 
above predicts that this does not become noticeable until temperatures below those of the clusters in their 
sample. 
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progressively smaller baryon fractions within the virial radius. The heating required to 
produce the observed entropy floor can be reduced to ~ 2.5 keV per particle for v = 0.25 
(implying a baryon fraction in 1 keV clusters nearly three times lower than in 10 keV clusters). 
It is implicitly assumed that the reduction in /viriai (by mass leaving the cluster, or simply 
expanding beyond rviriai) is a result of the heating: lower values of e can be accommodated 
if, instead, cooler clusters are assumed to be intrinsically gas-poor. However, the magnitude 
of the entropy gap between observed cool clusters and that expected assuming no heating 
and constant /viriai would require an extreme gas-fraction trend (PCN) of which there is no 
observational evidence (EF, Mohr et al. 1999). 

Finally, I have considered a family of mass sequences where the polytropic index in the 
heated model decreases by 50% for each factor of two decrease in mass {i.e. from n = 100 
to n ~ 11 as Tc - defined as the emission- measure-averaged ICM temperature within rsoo - 
varies from 10 to 1 keV), as might be appropriate if heating leads to more-nearly isentropic 
configurations (Balogh et al. 1999). Since as n decreases the ICM density distribution 
steepens, larger values of e - corresponding to ~ 4.5 keV per particle - are demanded by the 
observed entropy floor. 

One can conclude from the above that the observed entropy floor for cooler clusters is 
consistent with a mass sequence of polytropic equilibrium configurations heated by > 3 keV 
per particle at the cluster center (see Figure 2; 3 keV is about the temperature at which the 
observations start to deviate significantly from the predictions of models without heating). 
In the context of these models, this same amount of heat simultaneously accounts for the 
observed steepening of the relationship between ICM mass within rsoo and temperature 
(Figure 4). Note that the former is a local indicator of ICM expansion, while the latter is 
a global indicator. Also note that any contrivance meant to lower the required amount of 
energy injection (e.g., by reducing the intrinsic gas fractions of cooler clusters), will destroy 
the simultaneous match to the observed steepening of the relationship between total mass 
within rsoo and temperature (Figure 3) or, equivalently, the increasing departure from T^ = 
Tviriai implied by the observed a-T^ relation (Figure 5). 



4. Discussion 

Given the striking success of the e ~ 0.35 heated family mass sequence in reproducing the 
observed trends in global ICM properties - in particular their monotonically increasing de- 
parture with decreasing temperature from expectations of models without extra-gravitational 
heating - it is worthwhile to consider the implications for the relative gas distributions in 
hot and cool clusters, the connection between heating and metal enrichment, and the astro- 
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physical interpretation of these models in the context of hierarchical formation of clusters. 



4.1. Relative Gas Distributions 

Clearly, the existence of an entropy floor caused by heating implies that cooler clusters 
will have more extended ICM, i.e both lower central densities and shallower density slopes. 
And indeed this was discovered by PCN in their analysis of i? 05^4 T cluster data. The effect 
of heating on the gas distribution is perhaps most evident in plots of the cumulative gas 
fraction. Figure 6 shows a comparison of these for 1 and 10 keV clusters with and without 
heating (the latter from the simple heated family of Section 3.2). In the absence of heating, 
cooler clusters are expected to have significantly higher gas fractions at small radii (relative to 
the virial radius) because of their more concentrated mass distributions. Including heating, 
the opposite is very strongly the case: at v^qq (~ O.Srviriai), 1 keV clusters are about twice as 
gas-poor as 10 keV clusters (Figure 7), even though it is assumed that they have identical 
gas fractions at rvmai- A similarly strong trend is seen for the ICM density distribution slope, 
as shown by Figure 8 which plots Amage = 3dlnpgas/dlnr at rsoo versus T^ for mass sequences 
with and without heating. The density distribution is predicted to become flatter at all radii 
greater than O.lrvmai for cooler clusters. 

Whether these trends have been detected - and, if so, to what quantitative extent - 
is currently problematic. Since both the density slope and gas fraction vary (possibly non- 
monotonically) with radius, comparisons must be made at radii corresponding to identical 
density contrasts (e.g., rsoo). Thus the total mass must be estimated (for cooler systems this 
must done from the data itself, not from an extrapolation of the hot-cluster mass-temperature 
relation) in order to make a meaningful comparison of density slopes as well as gas fractions. 
Moreover, deriving the density distribution for the cooler clusters {i.e., groups) is greatly 
complicated by uncertainties in background subtraction, in treating the emission from indi- 
vidual galaxies, and possibly by departures from spherical symmetry ( pVIulchaey fc Zabludofl 



1998] , [Hwang et al. 19"99D . Analysis along the lines of EF (who restrict themselves to the 



more luminous, and therefore hotter, systems) on the PCN ensemble of surface brightness 
profiles could prove illuminating. Interestingly, extended elliptical galaxy X-ray halos, with 
kT^ ~ 0.7-1 keV, tend to have Amagc ~ 0.5, in accord with the heating sequence of Figure 8 
and their likely physical continuity with ICM ([Mathews fc Brighenti 1998|). 
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4.2. The Connection Between Heating and Metal Enrichment 

The presently most reliable diagnostic of supernova heating in clusters - that presum- 
ably occurred during star formation driven galactic outflows at early epochs - is the ICM 
silicon abundance. Only iron abundances are determined more accurately; however, the 
amount of energy ejection associated with the mass of Fe in the ICM is sensitive to the 
poorly determined (e.g., [Gibson, Loewenstein, fc Mushotzky 1997|) relative fractions of Fe 



originating from Type la (SNla) and Type 11 (SNll) supernovae. Not only have accurate Si 
abundances been determined for a large sample of clusters QFukazawa et al. 1999|) , but the 



explosion energy per Si yield is similar in standard SNla and SNll models (e.g., Gibson et 
al. 1997) so that the energy injection associated with the observed amount of Si is robust 
to uncertainties in the SNla/SNll ratio. The equivalent temperature of the injected energy 
can be expressed as 

fcT/^ ~ l.QZsiCsN keV, (32) 

where Zsi is the ICM Si abundance, relative to solar, and csn is the average injected energy 
per supernova (perhaps boosted by stellar winds from SNll progenitors; [Leitherer fc Heckman 
1995| ) in units of 10^^ erg. 



Since Zsi typically is ~ 0.7 for rich clusters QFukazawa et al. 1999|) , supernova energy 



injection would seem to fall short by a factor of ~ 3 of providing the required heating (this 
problem is exacerbated for groups where the Si abundance is often estimated to be lower 
by an additional factor of two; [Davis, Mulchaey, fc Mushotzky 1999| ; [Hwang et al. 19"99| ). 



However, one must keep in mind that the requirement of 3 keV per particle of heating is 
for the cluster center only. If this heating were uniformly present cluster-wide, the resulting 
expansion of the ICM would lead to an extremely steep decrease of /viriai with mass. The 
fact that this is not observed implies that the energy injection decreases outward from the 
center. The central concentration of SN energy injection may be reflected in abundance 
gradients that have now been measured with ASCA (e.g., [Ezawa et al. 1997[) . Additional 



heating at the cluster center could result from early AGN activity, although the mechanism 
for efficiently dissipating AGN energy flux in the ICM is unclear. 

The average cluster-wide heating per particle can be estimated as 5CrkTq/2fimp ~ 
0.38kTg/fj,mp, where Cr is deflned in equation (22) and F !^ 0.8. Thus the global-average 
heating is on the order 1 keV per particle. This is more in line with that expected from 
nucleosynthetic considerations, but requires a high energy deposition efficiency into the ICM 
or esN >> 1. The latter is possible if the initial mass function of SNll progenitors is skewed 
toward high masses (fTimmes, Woosley, fc Weaver 1996|) . This amount of energy injection 
exceeds the escape temperature of the lowest mass systems, and could lead to removal of 
some fraction of the ICM - depending on the detailed history of the heating with respect to 
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the dynamical evolution of the cluster. As discussed in the previous subsection, it is unclear 
whether lower mass systems have systematically lower gas fractions. 



4.3. Physical Interpretation and the Epoch of Heating 

The amount of heat characterizing the e ~ 0.35 heated family mass sequence corresponds 
to that required to transform the ICM from a non-heated equilibrium state, and is thus an 
upper limit since the gas may be pre-heated. If the heating occurs at sufficiently low density 
{i.e., prior to cluster collapse), the energetic requirements for obtaining the observed entropy 
minimum are greatly reduced. 

PCN argue that this must be the case. However, this is based in part on consideration 
of a SN rate (for the particular case of the Coma Cluster) derived from the present stellar 
content and the assumption of a standard initial mass function (IMF) that is well-established 
to fall far short of providing the observed heavy metal enrichment (e.g., [Loewenstein fc 



[Mushotzky 1996] ). Using standard SN yields, the Si enrichment, in solar units, can be 



expressed as 

Zs. = 0.2 C-^^) ^, (33) 

V MicM / 0.01' ^ ' 

where Mgtars and Micm are the total cluster stellar and ICM masses, respectively, and rjsN is 

the total number of SN per solar mass of present-day stars. PCN, assuming tjsn = 0.01 as is 

appropriate for a standard IMF, thus derive an equivalent supernova temperature ~ 5 times 

lower than the value in equation (32) that is based on the actual observed rich cluster Si 

abundance. The relative epochs of cluster assembly and ICM heating cannot be constrained 

solely on energetic grounds. 

In pre-heating models the ICM is heated to some initial temperature kT^ and high 
entropy state. The ICM in the lowest mass systems is spared any further entropy increase 
via accretion shocks but, instead, is accreted and compressed adiabatically - thus preserving 
the initial entropy at all radii. An entropy floor of the kind reported by PCN is predicted 
at low ICM temperature, and can be expressed as 

nZ^/^kT = lOOSo keV cm^ (34) 



e 



with So ~ 1. 

Consider two extreme variations on the pre-heating scenario. In the flrst (e.g., PCN), 
the universe as a whole is pre-heated to kT^, at a redshift of z^, and, therefore, an electron 
density of 

rie. = fibaryon^^^^^(l + Z,f = 1.2 10-^(1 + Z,)^ Cm'^ (35) 



-16- 



( [Ibzzi fc Norman 1999|) , where Pcrit,o is the critical density at zero redshift, and fibaryon = 
0.0255/i7o~^ as imphed by standard Big Bang nucleosynthesis. That is, 

kT^ = 2.45 10"-'^S„(1 + z^f keV. (36) 

Reproducing the observed luminosity- and mass-temperature relations requires /cT* ~ 0.5 
keV (Cavaliere et al. 1999), i. e. z^, k, 13 for So = 1. In fact, however, to obtain the change 
in slope of the L^-Tx relation (that is the transition to adiabatic compression) at the correct 
temperature requires S^ = 2.5-10 (|Tozzi fc Norman 1999|) , implying z* = 3.5-8. That is, 
Eq = 1 can be reconciled with the mass-temperature relationship only if the universe was 
heated to 0.5 keV at a redshift greater than 10 - and, in this case, the transition from shock 
heating to adiabatic compression occurs at too low a mass to be consistent with the L^-T^ 
relation. The transition occurs at the appropriate mass if the level of the entropy floor has 
been underestimated by a significant factor, however pre-heating redshifts greater than 3 
are still required if the mass-temperature relation is to be recovered as well. Moreover, such 
a high initial entropy is predicted to lead to baryon fractions in groups only ~ 0.2 of the 
universal value (Cavaliere et al. 1999), which could be problematic since there are group 
baryon fractions observed to exceed 10% ( pVIulchaey et al. 1998| , [Hwang et al. 1999| ). 



In the scenario of Balogh et al., each proto-cluster is heated at turnaround; i.e. they 
assume 

rie. = fibaryon^^^ ~ 8.3 10-^g{z,){l + z,f cm'^ (37) 

for Qo = 0.3 and h-ro = 1, where (7(2;*) = 1, 0.88, 0.81, and 0.75 at 2;* = 1, 2, 3, and 5, 
respectively. Thus, 

kT, ^ 1.2 lO-^g{z,y/^J:oil + z,f keV; (38) 

so that kT^ = 0.5 keV implies z* = 6.2. The models of Balogh et al. (1999) also require 
So > 1 (So ~ 3.7) to explain the low-temperature steepening of the L^-T^ relation. 

Both pre-heating scenarios require an initial entropy several times greater than the 
observed entropy floor as estimated by PCN. These can perhaps be reconciled if one recalls 
that the PCN entropy floor is deflned using a global average temperature. Pre-heating implies 
that accretion shocks become weak and the ICM increasingly isentropic at low temperature 
(well below where the entropy- and luminosity-temperature relations change slope, at kT^ ~ 
3 keV). In the limit of pure adiabatic compression, significant temperature gradients are 
expected. As a result the local entropy at Vc may be larger than the PCN estimate. 

One can consider pre- heated mass sequences of models as defined in Section 2, by fixing 
the polytropic index n at 3/2. Such models are normalized by fixing two of the following 
three parameters: So, the mass-averaged temperature in units of the virial temperature 
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Tm, and the gas fraction within the virial radius /vmai- A general feature of these mod- 
els is that the ICM temperature at r^ is about 50% greater than the emission-averaged 
temperature, and the entropy at r^ correspondingly higher than that estimated using the 
emission-weighted average. Results for an isentropic mass sequence with So = 2 (chosen to 
approximately reproduce the PCN entropy floor) and r^ = 1.05 are shown by the dotted 
curves in Figures 3 and 4 (sequences with constant /viriai have similar properties at low T^). 
The mass-temperature relation parallels that of the unheated sequence, shifted over by ~ 2 
in temperature, so that the observed relation is well- matched at low Tx (Figure 3). The 
ICM mass-temperature relation is reproduced as well for kT^ < 1.5 keV (Figure 4), and 
only starts to deviate strongly beyond 2 keV (since the gas fraction exceeds 0.5) where the 
isentropic assumption should start to break down. 

The temperature profile in the isentropic models are relatively steep (Figure 9). This 
can be compared to the observed gas distribution for HCG 62 - the best measured group 
temperature profile in terms of accuracy and extent. The density profile is flat as predicted, 
rie ~ r~^ ( [Ponman fc Bertram 1993]) ; but, for an isentropic distribution the temperature 
is expected fall to below 0.5 keV beyond 500 kpc - and does not ( Ib'inoguenov fc Ponman 
119991) . 

An assessment of the pre-heating scenario can now be summarized. If heating occurs 
prior to collapse, clusters with temperatures sufficiently below the turnover of the luminosity- 
(or, equivalently, the ICM mass-) and entropy-temperature relations are expected to be isen- 
tropic and display a steep trend of decreasing baryon fraction with decreasing temperature. 
The present observational evidence does not support this, although the existing data is not 
of sufficient quality to be conclusive. Moreover, to effect these turnovers at the appropri- 
ate mass scale requires an initial entropy at least twice the level of the observed entropy 
floor. In order to simultaneously explain the turnover in the mass-temperature relation the 
proto-cluster gas must be pre-heated to kT^ > 0.3 keV at redshift 

\ ' 'baryonPcrit l^* J / 

where the pre-collapse overdensity (/ie'^p^e*)/(^baryonPcrit(^*)) must - by definition - lie 
between one and the overdensity at turnaround (~ 15). In the context of the thermal 
history of baryons in the universe ( |Cen fc Ostriker 1999| ), such a level of pre-heating seems 
unlikely at that high a redshift. In fact the (mass-averaged) entropy in the baryon content 
of the universe is a strongly decreasing function of redshift, the universe having been both 
much cooler and more dense in the past, and does not approach the level of the observed 
cluster entropy floor until relatively recently. 

Thus, I would conclude that pre-heating models in their present preliminary form cannot 
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simultaneously explain all of the heating-induced departures from self-similarity observed in 
clusters - although a definitive assessment must await more refined pre-heating models and 
their comparison with more extensive and complete observations of low-mass systems. On the 
other hand, the success of the heated families of equilibrium polytropes in simultaneously 
explaining the deviations from self-similarity observed in clusters (Figures 2-5) suggests 
that substantial heating - certainly in the cluster core - occurred after the cluster was, at 
least partially, in place. Perhaps some of the heating was triggered by merging activity 
associated with the latter stages of cluster formation. It must be noted, however, that the 
metal enrichment (or, at least, the Fe enrichment) that presumably accompanies the heating 
cannot have occurred much more recently than z = 1 based on the lack of evolution of cluster 
abundances ( |Mushotzky fc Loewenstein 1997| ), and the properties of elliptical galaxies - the 
sites of primordial star formation generally presumed responsible for ICM heating and metal 
enrichment - as a function of redshift ( |Broadhurst fc Bouwens 1999| , and references therein). 

A detailed consideration of these issues requires a fully time-dependent approach (e.g., 
Kay fc Bower 1999|) , and must reconcile the required amount of heating with the ob- 



served properties of cluster galaxies (Wu et al., in preparation) - which may be problematic 
( IBrighenti fc Mathews 1999| ). 



4.4. Cool Cluster Caveats 

The breaking, by energy injection, of gravity-only self-similarity by an increasing amount 
for cooler clusters raises the following caveats. (1) In deriving the cluster mass function from 
X-ray luminosity (L^) or temperature functions, care must be taken to avoid using single 
slope Lx-Tx or M-T^ relations extending all the way down from high T^. (2) The gas fraction 
is a fairly steep function of radius for cool clusters (Figure 6) and values measured out to, e.g, 
rsoo, may not accurately reflect the global value. Groups may appear gas-poor in comparison 
to rich clusters while, in fact, their ICM are simply more extended with respect to their non- 
baryonic components. 

If hierarchical clustering correctly describes the formation of large scale structure, cooler 
clusters will tend to form over a wider range of redshifts (Balogh et al. 1999); and presumably, 
there will be greater variation in the relative timing of merging/collapse and star formation 
{i.e. cluster construction and heating). The large scatter observed in metal abundance and 
gas fraction (e.g., Davis et al. 1999) may be a reflection of this diversity. 



-19- 



5. Summary and Conclusions 

In this paper I liave attempted to construct a simple framework for understanding 
the relationship between the baryonic and non-baryonic components in galaxy clusters as 
manifest in the observed scaling of total mass, and ICM mass and central entropy with ICM 
temperature T^. 

Observed, as well as simulated, ICM density distributions of the most massive clus- 
ters are well-characterized by nearly isothermal hydrostatic configurations, with /viriai = 
0.155/i7o~^" and T^ = l.lTviriai, in an NFW potential. However, there are clear and system- 
atic departures from the expectations derived from simply extending application of these 
baseline models to cooler, less massive systems. These departures are in the sense that the 
ICM in less massive systems have an increasingly large excess of thermal energy (from the 
steepening of the total mass-temperature relation), and become increasingly less concen- 
trated - both globally (from the steepening of the total mass-temperature relation) and, 
more prominently, locally in the cluster core (from the flattening of the central entropy- 
temperature relation). 

This implies that substantial heating has occurred, heating with an increasingly pro- 
found effect on less massive systems. I considered transformations of the baseline models 
to families of new configurations that have specified amounts of additional heating at the 
cluster center (and possibly different gas fractions and polytropic indices), and found that 
mass sequences heated by a universal amount per particle can simultaneously reproduce the 
observed total mass-temperature, ICM mass-temperature, and central entropy-temperature 
trends to excellent precision. The required central heat input is substantial, > 3 keV per 
particle; however, the cluster-wide average value is ~ 1 keV per particle: just consistent with 
robust global nucleosynthetic constraints derived from rich cluster Si abundance measure- 
ments and standard supernova energies. The concentration of the required heat injection 
is reflected in the fact that the local (central) entropy-temperature relation shows a more 
extreme departure from the gravity-only prediction than the global ICM mass-temperature 
relation. 

The 3 keV per particle of heating required at the centers of clusters is an upper limit 
since it is derived assuming no pre-heating. This limit is derived from static considerations 
without recourse to a full dynamical treatment (with its many accompanying parameters), 
and becomes increasingly accurate the later the epoch of heating with respect to cluster 
virialization. PCN argue that the heating must occur prior to the epoch of cluster assembly. 
However, this is based in part on an underestimate of the available energy from supernova 
explosions by a factor of 3-5. If the universe as a whole is pre-heated, the total mass- 
temperature relation implies heating the baryons to ~ 0.5 keV at 2; > 3 which is unlikely given 
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that few objects with that virial temperature will have collapsed and most star formation 
is yet to occur. Pre-heating becomes more feasible the higher the overdensity of the proto- 
cluster at that epoch; however, pure pre-heating models generally require a higher entropy 
floor than is observed in order to reproduce the luminosity-temperature relation (equivalent 
to the ICM mass-temperature relation successfully reproduced by the heated mass sequence 
presented in Section 3.2). A strong prediction of pre-heating models is that the lowest mass 
clusters should be isentropic, a prediction difficult to assess with presently existing data 
although the entropy in the ~ 1 keV group HCG 62 appears to increase with radius. 

The form of the observed departures from self-similarity would seem to imply the ma- 
jority of ICM heating, i.e much of the massive star formation, was contemporaneous with 
- and perhaps caused by - the assembly of the cluster, occurring at redshifts between 10 
and 1. Comparison of more detailed evolutionary models that incorporate time- and space- 
dependent heating spread out over this epoch with improved X-ray observations - particu- 
larly of low-mass systems - will greatly illuminate the early star formation history of galaxies 
and the effect this has on the surrounding environment. 

I am grateful to U. Hwang for computing and providing the data points and errorbars 
in Figure 2, to D. Horner for providing the data points in Figures 3 and 5, and to A. Fabian, 
R. Mushotzky, P. Tozzi, and the anonymous referee for providing useful feedback. 
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Fig. 1. — Comparison of the density profile in the following hydrostatic polytropic models 
with simulated cluster /5-model fit (solid curve): n = 5 and r = 1.1 (dotted curve), n = 100 
and r = 1 (short-dashed curve), n = 100 and r = 1.1 (long-dashed curve), where n is the 
polytropic index and r the ratio of the emission-measure-weighted average temperature to 
virial temperature. 
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Fig. 2. — Entropy-temperature curves for heated models with n = 100, constant fvMai, 
and e = (solid curve) or e = 0.35 (dashed curve), where the entropy is defined as 
logio(/i;T^^e^^^'^)) "^e is in cm~^ and evaluated at r^, kT^ is the average (within rsoo) ICM 
temperature in keV, and h-j^ = 1. The filled triangle represents the value from numerical 
simulations without heating; the data points, representing averages in ten temperature bins, 
and errorbars, representing the dispersion within each bin, have been derived from ASCA 
data and kindly provided by U. Hwang. 
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Fig. 3. — Mass-temperature curves for e = and e = 0.35 models as in Figure 2, best-fit 
observed correlation from EF (stars), and isothermal /?-model masses from Horner et al. 
(1999; data points with errorbars), where the mass is the total gravitational mass integrated 
out to rsoo- Also shown is an isentropic mass sequence with mass- averaged temperature 
approximately equal to the virial temperature and kTrji^r'^l^ = 200 keV cm^ (dotted curve; 
see Section 4.3). 
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Fig. 4. — Same as Figure 3 for the ICM mass; best-fit observed correlation is from Mohr et 
al. (1999). 



-28- 



3.2 



o 



^ 



2.8 



2.6 



2.4 



1 — \ — I — I — I — I — I — I — I — I — I — r 



T — I — I — I — r 




0.2 0.4 0.6 0.8 

log <kT> (keV) 



Fig. 5. — Same as Figure 3 for (/cTviriai/yU^p)^''^; observed trend is best-fit a-T^ correlation 
from Bird et al. (1995); data points with errorbars are from Girardi et al. (1998). 
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Fig. 6. — Gas fraction profiles for heated models with n = 100, constant /viriai, and e = 
(solid curve: 10 keV cluster, dotted curve: 1 keV cluster) or e = 0.35 (dashed curve: 10 keV 
cluster, dot-dashed curve: 1 keV cluster). 
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Fig. 7. — Gas fraction at rsoo for heated models with n = 100, constant /vmab and e = 
(solid curve) or e = 0.35 (dashed curve). 
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Fig. 8. — Same as Figure 7 for /^image (three times the ICM density slope) at rsoo- 
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Fig. 9. — Mviriai = 2.3 IO^^Mq (Tvinai ~ 0.57 keV) cluster temperature profiles for e = 
(solid curve) and e = 0.35 (dashed curve) heated models featured in Figures 2-8, as well for 
the kT^rie'"^/^ = 200 keV cm^ isentropic model (dotted curve). The latter two models have 
emission-averaged temperatures ~ 1 keV, although the mass-averaged temperature in the 
isentropic model is, by construction, equal to 1.05Tviriai- 



